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The magnetic phase diagram of the perovskite-type Ti oxides as a function of the GdFeOs-type 
distortion is examined by using the Hartree-Fock analysis of a multiband d-p Hamiltonian from 
a viewpoint of competitions of the spin-orbit interaction, the Jahn- Teller (JT) level-splitting and 
spin-orbital superexchange interactions. Near the antiferromagnetic (AFM)-to-ferromagnetic 
(FM) phase boundary, A-type AFM [AFM(A)] and FM states accompanied by a certain type 
of orbital ordering are lowered in energy at large JT distortion, which is in agreement with 
the previous strong coupling study. With increasing the GdFe03-type distortion, their phase 
transition occurs. Through this magnetic phase transition, the orbital state hardly changes, 
which induces nearly continuous change in the spin coupling along the c-axis from negative to 
positive. The resultant strong two-dimensionality in the spin coupling near the phase boundary 
is attributed to the strong suppression of Tn and To, which is experimentally observed. On the 
other hand, at small GdFeC>3-type without JT distortions, which correspond to LaTiOs, the most 
stable solution is not G-type AFM [AFM(G)] but FM. Although the spin-orbit interaction has 
been considered to be relevant at the small or no JT distortion of LaTiOs in the literature, our 
analysis indicates that the spin-orbit interaction is irrelevant to the AFM(G) state in LaTiOs 
and superexchange-type interaction dominates. On the basis of further investigations on the 
nature of this FM state and other solutions, this discrepancy is discussed in detail. 

KEYWORDS: perovskite-type Ti oxides, GdFe03-type distortion, d-level degeneracy, d-type Jahn-Teller distortion, 
spin-orbit interaction, multiband d-p model 



§1. Introduction 

In transition-metal oxides, strong electron correlations 
often localize the 3d electrons and the system becomes 
an insulator (a Mott insulator). 1 ) These compounds 
have recently attracted considerable interest since they 
show rich magnetic and orbital phases. In particular, 
perovskite-type oxides RMO3, where R denotes a triva- 
lent rare-earth ion (i.e., La, Pr, Nd, Y) and M is a 
transition- metal ion (i.e., Ti, V, Ni, Cu) exhibit a va- 
riety of magnetic and electronic properties caused by an 
interplay of charge, spin and orbital degrees of freedom. 

The perovskite-type Ti oxide i?TiC>3 is a prototypical 
example. In these compounds, Ti 3+ has a t\ g configu- 
ration, and one of the threefold i2 ff -orbitals is occupied 
at each transition-metal site. They have attracted much 
interest since these systems show various magnetic and 
orbital ordcrings owing to the threefold degeneracy of 
the tig orbitals. It requires to take both spin and orbital 
fluctuations into consideration to explain competitions 
of such rich phases. Moreover, the spin-orbit interaction 
would make the magnetic and orbital structures more 
complicated since the t^g orbitals are strongly affected 
by the interaction. 

The crystal structure of i?TiC>6 is an orthorhombically 
distorted cubic-perovskite (GdFeC>3-type distortion) in 
which the TiC>6 octahedra forming the perovskite lat- 
tice tilt alternatingly as shown in Fig. 1. The magni- 
tude of the distortion depends on the ionic radii of the 
R ions. With a small ionic radius of the R ion, the 



lattice structure is more distorted and the bond angle 
is more significantly decreased from 180°. In LaTiC>3, 
the bond angle is 157° (a6-plane) and 156° (c-axis), but 
144° (afr-plane) and 140° (c-axis) in YTi0 3 2 >. The dis- 
tortion can be controlled by the use of the solid-solution 
systems Lai-j,Yj,Ti03 or in -RT1O3, by varying the R 
ions. In particular, by varying the Y concentration in 
Lai- a Y a Ti03, we can control the bond angle almost con- 
tinuously from 157° (y = 0) to 140° (y = 1). 

In YTiC>3, a G?-type JT distortion has been observed 
in which the longer and shorter Ti-0 bond lengths are 
~2.08 A and -2.02 A, respectively. 3 ) In the d-type JT 
distortion, the xy and yz orbitals are stabilized at sites 
1 and 3, and the xy and zx orbitals are stabilized at 
sites 2 and 4. On the other hand, LaTiOs exhibits no 
detectable JT distortion. 

Recently, the magnetic phase diagrams have been 
studied as functions of the magnitude of GdFe03-type 
distortion. 4_7 - ) In La-rich (y < 0.6) systems or in the 
compounds with large R ions, in which the GdFe03-type 
distortion is relatively small, an AFM ground state is re- 
alized. In particular, LaTiOs exhibits a AFM(G) ground 
state with magnetic moment of 0.45 /iR, which is strongly 
reduced from spin-only moment, and the Neel tempera- 
ture (Tn) is about 130 K. With increasing the Y concen- 
tration or varying the R site with smaller-sized ions (an 
increase of the GdFe03-type distortion), Tn decreases 
rapidly and is suppressed to almost zero, subsequently a 
FM ordering appears. In Y-rich systems and in YTi03 
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Fig. 1. GdFe03-type distortion. 



in which the GdFe03-type distortion is relatively large, 
the system shows a FM ground state. 

In order to elucidate these phase diagrams, model 
Hartree-Fock studies have been done previously. 8 ' 9 ^ In 
these weak coupling studies, it is claimed that in LaTi03 
with small GdFeC>3-type distortion, a AFM(G) state 
with the spin-orbit ground state is realized for the small 
or no JT distortion, and resultant unquenched orbital 
moment is considered to be consistent with the strong 
reduction of the moment. On the other hand, a FM state 
accompanied by an orbital ordering is realized in YTi0 3 
with large JT distortion. However, in these studies, the 
nature of the phase diagrams has not been elucidated 
sufficiently in the following sense. At first sight, we can 
expect the first-order transition between completely dif- 
ferent symmetry breaking in which Tn and Tc remain 
nonzero at the AFM-FM phase boundary. However, in 
the magnetic phase diagrams, Tn and Tc are strongly 
suppressed around the phase boundary. This strong sup- 
pression implies a continuous-type transition at T = 
and contradicts our naive expectation. This second-order 
like phase transition is not explained in these studies, and 
has been an issue of interest. 

Recently, in order to clarify this problem, effective 
Hamiltonian in the insulating limit has been applied 
to this system. 10 ' 11 ) According to these strong coupling 
studies, in the AFM phase near the AFM-FM phase 
boundary, an AFM(A) ground state is realized. This 
AFM (A) phase has not been studied in the previous weak 
coupling approach. Moreover, since the orbital state is 
strongly stabilized and changes only little through the 
transition, strong two-dimensionality in spin-coupling is 
predicted near the phase boundary and the strong sup- 
pressions of Tn and Tc are naturally understood. In 
these studies, a large JT distortion is assumed in or- 
der to focus on the situation near the AFM-FM phase 
boundary, and the spin-orbit interaction is neglected on 
the basis of the large energy-splitting due to the JT dis- 
tortion. In addition, the AFM(G) state in LaTi03 with 
small or no JT distortion has not been reproduced. 

However, the spin-orbit interaction may become rel- 
evant in the systems with small or no JT distortion 
such as LaTi03. While the spin-orbit interaction can 
be neglected if the JT distortion is large, we can expect 
a strong competition between the spin-orbit interaction 



and JT level-splitting with decreasing the JT distortion. 
With sufficiently small JT distortion, the system may 
well be described by the spin-orbit ground state. Be- 
sides, successive spin-orbital superexchange interactions 
may dominate over the spin-orbit interaction even with- 
out JT distortion. At this stage, it is an issue of impor- 
tance to examine the phase diagrams from a viewpoint 
of their competitions. In the weak coupling approach 
in which transfers of electrons and spin-orbit interaction 
are treated in a non-perturbative manner, we can con- 
sider both effects on an equal footing. This approach 
is appropriate for a systematic study on the interplay of 
them. In these senses, the weak-coupling and the strong- 
coupling studies are complementary to each other, and 
analysis from the weak coupling approach is important. 

In this paper, we investigate the magnetic phase dia- 
grams by using the Hartree-Fock analysis of the multi- 
band d-p Hamiltonian. We study the magnetic and or- 
bital states as functions of the GdFeC>3-type and d-type 
JT distortions. Since effects of both electron trans- 
fers and spin-orbit interaction are taken into account 
on an equal footing, this model is appropriate for a 
study on the competitions of spin-orbit interaction, JT 
level-splitting and superexchange interactions. The weak 
coupling treatment does not properly reproduce the en- 
ergy scale of the superexchange interaction J defined in 
the strong coupling region, where J is proportional to 
t 2 /U with t and U being typical transfer and on-site 
Coulomb repulsion. However, the physics contained in 
the reproduction of the superexchange interaction with 
AFM and/or antifcrro-orbital (AF-orbital) is expected 
to be adiabatically connected with the SDW type sym- 
metry breaking in the weak-coupling Hartree-Fock so- 
lution. Therefore, we will refer the stabilization of the 
SDW (or orbital density wave) type solution with AFM 
(or AF-orbital) symmetry breaking to the superexchange 
mechanism. 

Pioneering works by using this method have already 
been done by Mizokawa and Fujimori. 8 ' 9 ) However, con- 
cerning the region of small GdFe03-type distortion, we 
have come to a different conclusion by studying orbital- 
spin states which they have overlooked. We show that 
in the small GdFe0 3 -type distortion without JT level- 
splitting, a FM spin state accompanied by an AF-orbital 
ordering is stabilized by the energy gains of both spin- 
orbit and superexchange interactions. In this FM solu- 
tion, the spin-orbit ground state is not realized at certain 
sites, which suggests that the spin-orbital superexchange 
interactions due to the electron transfers dominate over 
the spin-orbit interaction even without JT distortion. 
In the previous studies, AFM(G) state with spin-orbit 
ground state has been claimed to be stabilized without 
JT distortion. However, in these studies, the stabiliza- 
tion of this AFM(G) state is concluded only from com- 
parison of the energies between this AFM(G) solution 
and a FM solution with higher energy, and our FM solu- 
tion is ignored. Our FM state has not been reproduced 
so far, and is studied for the first time by our weak cou- 
pling approach. We conclude that the AFM(G) state 
in LaTi03 does not accompany with spin-orbit ground 
state, and there exists another origin for its emergence. 



Spin and Orbital States and Their Phase Transitions of the Perovskite-Type Ti Oxides: Weak coupling Approach 



Recent neutron-scattering experiment shows the spin- 
wave spectrum of LaTiC>3 well described by a spin- 1/2 
isotropic Heisenberg model on the cubic lattice, and ab- 
sence of unquenched orbital momentum. 12 ' 1 This result 
also seems to contradict the naive prediction of spin- 
orbit ground state with no JT distortion. By studying 
a model including the spin-orbit interaction, we propose 
some statements on this experimental result. 

Moreover, we apply this method to the systems near 
the AFM-FM phase boundary for the first time. We 
show that the strange behavior of the magnetic phase 
transition is well described on the basis of JT ground 
state when we consider the experimentally observed large 
JT distortion. The results on the properties and na- 
ture of the phase transition are in agreement with those 
obtained by the previous strong coupling approaches, 
which indicates its validity irrespective of the coupling 
strength. 10 ' n ) In addition, we study a magnetic and or- 
bital phase diagram in the plane of the GdFeC>3-type and 
c?-type JT distortions in order to examine how extent the 
physics of AFM-FM phase transition in strong coupling 
limit survives when the JT level-splitting competes with 
the spin-orbit interaction. 

The organization of this paper is as follows. In § 2, we 
introduce the multiband d-p Hamiltonian to describe the 
realistic systems of the perovskite-type Ti oxides. In § 3, 
numerical results calculated by applying the unrestricted 
Hartree-Fock approximation are presented. Section. 4 is 
devoted to the summary and conclusions. 

§2. Multiband d-p model 

We employ the following Hamiltonian: 

H dp = H da + H p + H t dp + H tpp + H h + i?on-site, (1) 

with 

H d Q = ^ £ ddi.i 1<7 d a,z 7 cr, (2) 
H p = £ P Pl,jlaP«,jl<r, (3) 

H td p= Cy,a'^4,i 7CT P«',j;<T+h.C., (4) 

a,i,"y,ot f 1 j,l,(T 

H tpp = 1 %jl,a'j'l'PLjl*Pa'j'l'<r + h.C., (5) 

aj.l.a' ,j' ,1' ,a 

Hh ^ ' hj t7 ^i rT id a ^ cr d ctt i-yi tT i, (6) 

-Hon-sitc = H u + H u i + Hj + Hj' , (7) 

where d^ ■ is a creation operator of an electron with 
spin a(= j, J.) in the 3d orbital 7 at Ti site i in the a-th 
unit cell, and p^ a . J<T is a creation operator of an electron 
with spin a(= |, J.) in the 2p orbital I at oxygen site j in 
the a-th unit cell. Here, Hdo and H p stand for the bare 
level energies of Ti 3d and O 2p orbitals, respectively. 
H t dp and H tpp are d-p and p-p hybridization terms, re- 
spectively. Hh denotes the crystal field and spin-orbit in- 
teraction represented by the parameter £ = 0.018 eV. 13 ^ 
The term fl" on _ s i te represents on-site d-d Coulomb in- 
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teractions. t ap llaljl and t pp - lalj , v are nearest-neighbor 
d-p and p-p transfers, respectively, which are given in 
terms of Slater-Koster parameters V v dn, V p da, Vppir and 
Vppcr. 14 * 1 i? on _ s ite term consists of the following four con- 
tributions: 

Hu = ^ U ^a.i 1 ]^a.i^\d i a il ^d a , il i 1 (8) 
a, i. 7 

Hu' = u d) a il<y d a ^ ia d} a i y a ./d oli ij' cr i , (9) 

a,«,7>7 / ,cr,(j / 

Hj = jd^^^da^ia-d^^/^/da^ijtj/, (10) 

a,2,7>7 / <7,<r' 

Hj' = Y j'dl^da^'^dl^^a.iYi, (11) 

a, i, 7^7' 

where H u and H u i are the intra- and inter-orbital 
Coulomb interactions and Hj and Hji denote the ex- 
change interactions. The term Hj is the origin of the 
Hund's rule coupling which strongly favors the spin align- 
ment in the same direction on the same atoms. These 
interactions are expressed by Kanamori parameters, u, 
u', j and j' which satisfy the following relations: 15 ' 16 - ) 

20 

u = U+-j, (12) 

u' = u- 2j, (13) 

3=3'- (14) 

Here, U gives the magnitude of the multiplet-averaged 
d-d Coulomb interaction. The charge-transfer energy A, 
which describes the energy difference between occupied 
O 2p and unoccupied Ti 3d orbitals, is defined by U and 
energies of the bare Ti 3d and O 2p orbitals e d and e p as 
follows, 

A = e d + U-e p . (15) 

The values of A, U and V v da are estimated by the cluster- 
model analyses of valence-band and transition-metal 2p 
core- level photoemission spectra. 17 ' 18 ) We take the val- 
ues of these parameters as A = 7.0 eV, U = 4.0 eV, 
Vpda — — 2.2 eV and j = 0.64 eV throughout the present 
calculation. The ratio Vpda/Vpdn is fixed at —2.18, 
and Vp pa and V pp7T at 0.60 eV and —0.15 eV, respec- 
tively 19 " 21 ). The effects of the GdFe0 3 -type distortion 
and the d-type JT distortion are reflected on the hopping 
integrals. The GdFe0 3 -type structure is orthorhombic 
with orthogonal a-, b- and c-axes which can be obtained 
by rotating the four octahedra in the unit cell. Let us 
represent the four octahedra in the unit cell as site 1 , site 
2, site 3 and site 4 as shown in Fig. 1. Here, we simulate 
the GdFe03-type structure by tilting the TiC>6 octahe- 
dra by +8 and —6 about the (1,1,1) and (—1, —1, 1) axes 
with respect to the x,y and z axes. The magnitude of 
the GdFe03-type distortion is expressed by the bond an- 
gle. The magnitude of the JT distortion can be denoted 
by the ratio [V^/V^J 1 / 3 ; here, V p s da and are the 
transfer integrals for the shorter and longer Ti-0 bonds, 
respectively. The value for YTiOs estimated by using 
Harrison's rule takes ^1.040 19 '. This large JT distortion 
is also considered to be realized near the AFM-FM phase 
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Fig. 2. Energies of various spin and orbital configurations relative 
to that of AFM(A) state as functions of the Ti-O-Ti bond angle 
in the case of [V° d JV' d jV 3 = 1.040. 



boundary. 

We can rewrite the Hamiltonian in the k-space form 
by using the following Bloch-electron operators, 



4. 



1 yjk-R c 



(16) 



(17) 



P k,jla ^ N 

v a 

where k labels the wave vector in the first Brillouin zone. 

§3. Results and Discussions 

In this section, we present the numerical results cal- 
culated by applying the unrestricted Hartree-Fock ap- 
proximation to the multiband d-p model introduced in 
the previous section. In our calculations, we have con- 
centrated on uniform solutions. At this stage, the order 
parameters can be written as, 



(4 



k 



N 



(18) 



which are to be determined self-consistently. We have 
taken 512k points in the first Brillouin zone of the 
GdFeC>3-type structure and iterated the self-consistency 
cycle until the convergence of all the order parameters 
within errors of 1 x 10~ 4 . It should be noted that the 
basis of the Ti 3d orbitals are defined by using x-, y-, and 
z-axes attached to each TiOg octahedron in this paper. 

First, in order to focus on the situation near the AFM- 
FM phase boundary, the magnitude of the JT distortion: 
Wpda/VpdaY^ ls nx ed at 1.040, which is considered to be 
realized around the AFM-FM phase boundary. 

In Fig. 2, relative energies of various spin and orbital 
configurations are plotted as functions of the Ti-O-Ti 
bond angle from 157° to 140°. In the small GdFe0 3 - 
type distortion, a FM solution with (yz, xy, xy, zx)-type 
orbital ordering in which site 1,2,3 and 4 are dominantly 
occupied by yz, xy, xy and zx, respectively (FM1 solu- 
tion) is stabilized (see Fig. 3) since the FM state with the 
orbital configuration in which the neighboring occupied- 
orbitals are approximately orthogonal (AF-orbital order- 
ing) is favored both by transfers and by the exchange 
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Fig. 3. The orbital occupation in the majority spin states of FM1 
state as a function of the Ti-O-Ti bond angle in the case of 
[V, IV 1 . 1 1/3 = 1.040. 
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interaction j. A AFM(G) solution with (yz,xy,xy, zx)- 
type orbital ordering [AFM(G)1] has much higher en- 
ergy. However, it should be noted that the present cal- 
culations are carried out in the case of large JT distor- 
tion so that the obtained FM1 solution with the small 
GdFe03-type distortion does not necessarily contradict 
the emergence of AFM(G)-ground state in LaTiOa with 
no JT distortion. 

As the GdFe03-type distortion increases, the 
(yz, xy, xy, zx)-type orbital state becomes unstable. In- 
stead, the solutions with the orbital state in which 
xy orbital is mixed into the occupied yz and zx or- 
bitals [(yz, zx, yz, zx)-type orbital state] become stable 
(see Fig. 4). By moderately increasing the distortion, 
AFM(A) state with (yz, zx,yz, zx)-type orbital order- 
ing is stabilized relative to FM1 solution. With fur- 
ther decreasing of the bond angle, the FM state with 
(yz, zx, yz, zir)-type orbital ordering (FM2 solution) is 
stabilized. The AFM(G) solution with (yz, zx, yz, zx)- 
type orbital ordering [AFM(G)2] has much higher en- 
ergy relative to the other solutions. The AFM(A) to 
FM2 phase transition occurs at ZTi-0-Ti~ 142°. These 
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Fig. 4. The orbital occupation in the majority spin states of (a) 
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Fig. 5. Magnetic moment of various spin and orbital states 
as functions of the Ti-O-Ti bond angle in the case of 
[V^/V^J 1 / 3 = 1.040. In the FM1 state, the magnitude of 
the magnetic moment is different between sites 1, 4 and sites 2, 
3. 



AFM(A) and FM2 states are expected to be realized in 
the systems which are located near the AFM-FM phase 
boundary. In addition, we note that in the large JT 
distortion of 1.040, the spin-orbit ground state does not 
have any stable solutions. 

In Fig. 5, we have plotted the magnetic moment of 
the various spin and orbital solutions as functions of the 



Ti-O-Ti bond angle. With this large JT distortion, the 
orbital angular momentum is mostly quenched and the 
magnetic moment basically consists of the spin-only mo- 
ment. This indicates that the effect of the spin-orbit in- 
teraction can be neglected and the system near the phase 
boundary is well described by the JT ground state. 

Under this large d-type JT distortion, the occupation 
of the higher ti g orbitals at sites 1, 2, 3 and 4 are close 
to zero so that the occupied orbitals at each site can 
be expressed by the linear combination of the twofold 
degenerate lowered orbitals, approximately. In addition, 

since the order of the indirect (i-d-transfers mediated by 

v 2 

the O 2p orbitals is — ^ ~ 0.2 eV and sufficiently small 
compared with U, the fe-dependence of the coefficients 
for the linear combinations can be neglected. So that, 
we can express the occupied orbitals at i-th site by the 
coefficients C a j 7<T as, 



7 



C a .i ia \l > 



(19) 



Here, ^Z 7 ' denotes the summation over the twofold low- 
ered orbitals in the d-type JT distortion at z-th site, 
namely, xy and yz orbitals at sites 1 and 3, and xy and 
zx orbitals at sites 2 and 4. Since the spin-orbit inter- 
action is not effective under the large JT distortion, the 
imaginary parts of the coefficients are negligible. At this 
stage, we define the absolute values of the coefficients in 
the normalized form as, 



\C, 



a,vy(T | 



(20) 



The orbital states realized in the AFM(A), AFM(G)2 
and FM2 solutions can be specified by using angles 

#AFM(A), #AFM(G)2 an d #FM as, 

site 1; cos9 x \xy > +si~n9 x \yz >, 
site 2; cos 9 x \xy >+ sin 9 x \zx >, 
- cos 9 x \xy > + sin.8 x \yz >, 



site 
site 



4; — cos 9 x \xy > + sin 9 x \zx >, (21) 



where x = AFM(A), AFM(G)2 and FM2. In Fig. 6, the 
angles for the AFM(A), AFM(G)2 and FM2 solutions 
are plotted as functions of the Ti-O-Ti bond angle. In 
AFM(A) state, the sites 1, 2, 3 and 4 are occupied by 
c\yz + c%xy, c\zx + c%xy, c\yz — c^xy and c\zx — C2xy 



1), respectively. In particular, the difference 



between c\ and C2 tends to be small with increasing the 
GdFe03-type distortion, and both c± and ci take ap- 
proximately 1/V2 for the large distortion. Moreover, 
the similar orbital state is also realized in FM2 state, 
and both c\ and C2 also take approximately l/v2- This 
orbital state is in agreement with previous theoretical 
predictions, 9, u ' 22, 23 - ) and is observed experimentally in 
YTi0 3 . 24_27) The difference between 9 



AFM(A) 



and 9 



FM2 



is very small, especially in the largely distorted region or 
near the AFM(A)-FM2 phase boundary. This indicates 
that the orbital ordering hardly changes through the 
magnetic phase transition. Then, the AFM(A)-to-FM2 
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Fig. 6. The orbital structures in AFM(A), AFM(G)2 and FM2 
solutions as functions of the Ti-O-Ti bond angle. The differ- 
ence between those of AFM(A) and FM2 is considerably small, 
particularly for the large GdFe03-type distortion or near the 
AFM(A)-FM2 phase boundary. 
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phase transition is identified as a nearly continuous one 
with a tiny jump in the spin-exchange interaction along 
the c-axis from positive to negative and it takes approxi- 
mately zero at the phase boundary. On the contrary, the 
FM spin-exchange interaction is constantly realized in 
the a6-plane. The resultant strong two-dimensionality in 
the spin coupling can cause the strong suppression of Tn 
and Tc near the phase boundary. These are all in agree- 
ment with the previous strong coupling studies 10, n ) and 
indicate that these results are valid even at a realistic 
and intermediate coupling strength. 

We can also specify the orbital state realized in the 
FM1 solution by using two angles d\ and 9 2 as follows, 

site 1; cos#i|yz > + sm0i\xy >, 

site 2; cos 9 2 \zx >+ sin 02\xy > , 

site 3; — cos 9 2 \yz > + sin9 2 \xy >, 

site 4; — cos 9\ \zx > + sin6>i \xy > . (22) 

In Fig. 7, the angles 9\ and 9 2 are plotted as func- 
tions of the Ti-O-Ti bond angle. In the small GdFe03- 



Fig. 8. Energies of various spin and orbital configurations relative 
to that of AFM(G)3 state as functions of the Ti-O-Ti bond angle 
in the case of [V^/V^J 1 / 3 = 1.000. 



type distortion (ZTi-0-Ti~ 157°), almost complete 
(yz, xy, xy, zx)-type occupation is realized. In this or- 
bital ordering, the neighboring occupied-orbitals are ap- 
proximately orthogonal and electron transfers from the 
occupied orbitals are restricted to neighboring unoccu- 
pied orbitals. However, with increasing the GdFe03-type 
distortion, the occupations of the xy, zx, yz and xy or- 
bitals gradually increase at sites 1, 2, 3 and 4, respec- 
tively (see also Fig. 3). Therefore, both 9\ and 9 2 tend 
to become close to 45°, and the orbital state in the FM1 
solution become similar to that in the FM2 solution as 
the GdFe03-type distortion increases. As a result, FM1 
solution in the large GdFe03-type distortion is similar 
to that of FM2 so that the energy difference between 
FM1 and FM2 solutions is small in the largely distorted 
region. This indicates that the (yz, zx,yz, zx)-type or- 
bital ordering realized in the AFM(A) and FM2 states is 
strongly stabilized for the large GdFe03-type distortion. 
In addition, we note that with the large JT distortion 
of 1.040, the spin-orbit ground state does not have any 
stable solutions. 

We next fix the magnitude of the JT distortion: 
i V pda/ V pdJ 1/3 at i- 000 ( Le - no JT distortion) in order 
to focus on the situation realized in LaTiOs. In Fig. 8, 
relative energies of various spin and orbital configura- 
tions are plotted as functions of the Ti-O-Ti bond an- 
gle. Without the JT distortion, AFM(G)1 and AFM(G)2 
states have no stable solutions. 

So far, in the small GdFe03-type distortion, the AFM 
state with spin-orbit ground state [AFM(G)3], out of 
which two states with antiparallel spin and orbital mo- 
ment, -^(yz + izx) | and -^(yz — izx) J. are alternating 
between nearest neighbors, is considered to be stabilized 
both by the spin-orbit interaction and by the superex- 
change interactions. However, though this AFM(G)3 
state is lower in energy relative to the FM3 state in which 
-^(yz + izx) T (^(yz ~ izx) |) is occupied at each site, 
a FM state with AF-orbital ordering (FM1) is always 
lower in energy as compared with AFM(G)3 and FM3 
solutions. This indicates that spin-orbital superexchange 
interactions caused by electron transfers dominate over 
the couplings of the spin and orbitals due to the spin- 
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Fig. 9. Magnetic moment of various spin and orbital states 
as functions of the Ti-O-Ti bond angle in the case of 

l V pd<r/ V pdJ 1/3 = L00 °- In the FM1 state > the magnitude of 
the magnetic moment is different between sites 1, 4 and sites 2, 
3. 



orbit interaction, and the spin-orbit interaction does not 
play a role in the emergence of AFM(G) state in LaTiC>3. 

In the FM1 state, the sites 1, 2, 3, and 4 are approx- 
imately occupied by ^(u z + izx) 1\ xy T> xy T an( l 
-^(yz + izx) T, respectively. At sites 1 and 4, the spin- 
orbit ground state with antiparallel spin and orbital mo- 
ment is realized. Since the neighboring -j^{yz + izx) f 
and xy are approximately orthogonal, AF-orbital order- 
ing accompanied by the spin-orbit ground state is real- 
ized in FM1 solution. Consequently, this FM1 state is 
strongly stabilized both by the spin-orbit interaction and 
by the spin-orbital superexchange interactions. 

In Fig. 9, the magnetic moment of various spin and 
orbital solutions are plotted as functions of the bond 
angle. In AFM(G)3 and FM3 with spin-orbit ground 
state, the magnetic moment is strongly reduced from the 
spin-only moment due to the antiparallel contribution of 
the unquenched orbital moment while those of AFM(A) 
and FM2 with JT ground state take approximately unre- 
duced values. In FM1, reduced ordered moment is real- 
ized at sites 1 and 4 with the spin-orbit ground state 
while the moments are not so reduced at sites 2 and 3 in 
which xy orbital is dominantly occupied. 

The strong stabilization of the FM1 state in the small 
GdFeOvtype distortion with no JT distortion indicates 
that the spin and orbital states in LaTiOa can not be de- 
scribed by the spin-orbit ground state. In addition, there 
exists a discrepancy between the calculated energy dif- 
ference and that expected from experimentally obtained 
Tn of ~130 K. We expect the energy difference between 



FM and AFM(G) solutions per unit cell from Tn in the 
following way. First, we can naively estimate the spin- 
exchange constant J in LaTiOa based on a comparison 
of Tn with the numerical study on the spin- 1/2 Heisen- 
berg model on a cubic lattice as J = /cb?n/0.946 ~ 
12 meV. 28 ) Then, a bond -energy difference between FM 
and AFM spin configurations per Ti-Ti bond is J/2, and 
there are 12 Ti-Ti bonds in the unit cell so that we can 
expect that AFM(G) solution is lower than FM solution 
in energy by 6J ~ 72meV within the Hartree-Fock ap- 
proximation. However, this value is considerably large as 
compared with the characteristic order of the calculated 
energy difference even if the spin-orbit ground state is 
realized in LaTi03. (For instance, the energy difference 
between AFM(G)3 and FM3 is ~ 1 meV per unit cell.) 
This discrepancy can not be explained within the error 
bars of the parameters estimated from the analyses of 
photocmission spectra so that we can conclude the spin- 
orbit interaction is irrelevant to the AFM(G) state in 
LaTi0 3 . 

Here, a question arises: why is the ordered moment 
reduced from 1 /ib so strongly if the spin-orbit inter- 
action can not be its origin? Recent optical measure- 
ment shows that LaTiOa has a considerably small opti- 
cal gap of ~ 0.1 eV in the vicinity of the metal-insulator 
(M-I) phase boundary with strong itinerant character. 6 ) 
Therefore, in this system, we expect that some amount 
of charge and spin fluctuations remain. The reduction 
of the magnetic moment may easily be attributed to this 
itinerant fluctuation. 29 ) For instance, in 2D case, the or- 
dered moment ~ 0.6 /^b for the Heisenberg model dimin- 
ishes to < 0.2 (Ub for U = 4 Hubbard model due to the 
itinerant fluctuation accompanied by the double occu- 
pancy. 29 ' This strong reduction of the ordered moment 
with charge fluctuations is also obtained for Hubbard 
model with next-nearest neighbor transfers in recent nu- 
merical study. 30 ' Within the Hartree-Fock approxima- 
tion, the ordered moment is equivalent to the local mo- 
ment so that the reduction of the moment can not be 
reproduced. However, we can expect this reduction irre- 
spective of the dimensionality in an insulator with small 
insulating gap near the M-I phase boundary. Conse- 
quently, though the spin moment within the spin-wave 
theory takes ~ 0.844 /xb and the reduction due to the 
quantum effects is small in 3D spin-wave approximation, 
the ordered moment would easily diminishes to ~ 0.45 
/xb in LaTiOs with the strong itinerant character and 
large expectation value of the double occupancy when 
charge fluctuations are properly taken into account. 

In Fig. 10, we show the magnetic and orbital phase 
diagram in the plane of the GdFe03-type and d-type 
JT distortions. In the region of [V^/V^] 1 / 3 > 1.027, 
AFM(A)-FM phase transition occurs as increasing the 
GdFe03-type distortion. In the small GdFe03-type dis- 
tortion, FM state with AF-orbital ordering is stabilized 
in the whole range of [V^/V^J 1 / 3 . In particular, in the 
small JT distortion, only FM1 state is stabilized. In FM1 
state with no JT distortion, sites 1, 2, 3 and 4 are occu- 
pied by -^(yz + izx) t, xy f, xy | and ^{yz + izx) |, 
respectively, and though the spin-orbit ground state is 
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Fig. 10. Magnetic and orbital phase diagram in the plane of the 
GdFe03-type and d-tjpe JT distortions. 



realized at sites 1 and 4, xy |-occupancy is favored at 
sites 2 and 3 by the spin-orbital superexchange interac- 
tions. In addition, AFM(G) phase does not exist even 
for the small JT distortion. 

§4. Summary and Conclusions 

In this paper, we have studied the magnetic and or- 
bital states and their phase transitions of the perovskite- 
type Ti oxides by using the multiband d-p Hamiltonian. 
In this Hamiltonian, effects of both electron transfers 
and spin-orbit interaction are considered. By apply- 
ing the unrestricted Hartree-Fock approximation to this 
Hamiltonian, we have investigated the orbital-spin states 
as functions of the magnitudes of GdFeC>3-type and d- 
type JT distortions from a viewpoint of competitions of 
the spin-orbit interaction, JT level-splitting and spin- 
orbital superexchange interactions. These competitions 
are characteristic in ti g systems such as titanates in con- 
trast with e g systems such as manganites since the spin- 
orbit interaction strongly affects the t 2g orbitals rela- 
tive to the e g orbitals and JT coupling is rather weak 
in t2 g systems while the coupling almost always domi- 
nates over the spin-orbit interaction in e g systems. Our 
model and approach which treat the electron transfers 
and the spin-orbit interaction on an equal footing and 
in a non-perturbative manner are appropriate for the 
study of the competitions. We expect that the physics of 
AFM or AF-orbital ordering with superexchange mech- 
anism in the strong-coupling region is connected adia- 
batically with the SDW-type symmetry breaking in the 
weak-coupling Hartree-Fock solutions. So that, we have 
referred the stabilization of the SDW (or orbital density 
wave) type solution with AFM (or AF-orbital) symmetry 
breaking to the superexchange mechanism. 

In the perovskite-type Ti oxides, the transfers of elec- 
trons on Ti 3d-orbitals are governed by supertransfcr 
processes mediated by the O 2p states. We can calcu- 
late the nearest-neighbor and next-nearest-neighbor d-d 
transfers (t and t', respectively) by using perturbational 
expansions with respect to d-p and p-p transfers which 
are determined by using Slater-Koster parameters. A 
tight-binding (TB) Hamiltonian with thus obtained t 



and t' well reproduces the band structure obtained in 
LDA calculations. 31 -* The characteristic perturbational 
processes for t and t' are mediated by one O ion and by 
two O ions, respectively. The order of t and t' are t^ d /A 
and tp d t pp /A 2 with t p d and t pp being characteristic d-p 
and p-p transfers, respectively. In these compounds, the 
order of t pp / A is about ~ 0.05 at most so that t' is much 
smaller than t. Actually, the band structure calculated 
by using TB model with both t and t' is almost the same 
as that obtained by using the model with only t, particu- 
larly in the i2trband dispersions. When t' is negligible as 
compared with t, we can expect that considerable degree 
of nesting remains. Consequently, in these system, the 
Hartree-Fock calculation can give reliable results for the 
AFM and AF-orbital type symmetry breaking. 

Similar weak coupling approach has already been ap- 
plied to both end compounds LaTiOa and YTi03 by Mi- 
zokawa and Fujimori. 8,9 ) On the other hand, in this pa- 
per, the systems located near the AFM-FM phase bound- 
ary are studied by this method for the first time. More- 
over, by studying a FM state with lower energy which 
they overlooked, we conclude that the spin-orbit interac- 
tion can not be an origin for the AFM(G) state in LaTiOa 
in contrast with their conclusion. The conclusions of this 
paper are as follows. 

In the region of large JT distortion, the spin-orbit in- 
teraction is dominated by the JT level-splitting and the 
system is well described by the JT ground state. In 
this region, the AFM(A)-to-FM phase transition occurs 
with increasing the GdFe03-type distortion. Through 
this phase transition, the orbital state changes negli- 
gibly in agreement with the previous strong coupling 
studies. 10,11 ^ The negligible change in the orbital state 
through this AFM(A)-FM phase transition causes a 
nearly continuous change in the spin-coupling along the 
c-axis, and we can attribute the strong suppressions of 
Tn and Tc to the resultant two-dimensionality in the 
spin coupling near the phase boundary. The orbital 
states obtained in the FM2 and AFM(A) solutions are 
in agreement with those obtained by the previous strong 
coupling approaches, 10, 11 - ) which indicates the validity of 
the results even at a realistic and intermediate coupling 
strength. Actually, this orbital state has already been 
observed in YTi03. 24 ~ 27 ) We expect that a similar or- 
bital ordering may be observed in the compounds near 
the AFM-FM boundary such as SmTi0 3 , GdTi0 3 and 
Lai-yYyTiOa (y ~ 0.3). Recent resonant x-ray scatter- 
ing study shows that the orbital states in SmTiOs and 
GdTiOa have twofold symmetry similarly to YTi03, and 
this seems to be in agreement with our result. 32 ) Here, 
we note that neutron scattering experiment reveals that 
the magnetic structure of the Ti sites in SmTiOs is not 
AFM(A) but AFM(G). 33 ) In addition, though SmTi0 3 
is located near the phase boundary, Tn of ~ 50 K is 
somewhat high relative to the previous theoretical pre- 
diction. 10, n ) In SmTi0 3 , there exist magnetic moments 
on the Sm sites, and Sm-Ti spin-coupling may be im- 
portant for its magnetic properties while our model does 
not take the orbital and spin degrees of freedom on the 
R sites into account. However, our model can well de- 
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scribe the orbital-spin states and their phase transitions 
of LaTi03, YT1O3 and Lai-^Y^TiOs systems with no 
magnetic moments on La and Y sites. Moreover, since 
the orbital state near the AFM-FM phase boundary is 
strongly stabilized irrespective of the spin structure as 
shown in both our weak-coupling and previous strong- 
coupling studies, the similar (yz, zx, yz, zx)-type orbital 
state is also expected to be realized in SmTi03 though 
the magnetic structure is AFM(G) due to the Sm-Ti 
spin-coupling. 

Without a JT distortion, owing to both spin-orbit 
and spin-orbital superexchange interactions a FM state 
with the spin-orbit ground state accompanied by an 
AF-orbital ordering [(^(u z + izx) f, xy f, xy |: 
-^(j/z + izx) t)-orbital ordering] is stabilized relative to 
the other solutions. This FM solution can not be repro- 
duced by the previous strong coupling approach in which 
the spin-orbit interaction is neglected in the large JT 
distortion, and is studied for the first time by our weak 
coupling approach. In addition, AFM(G) state is higher 
in energy and has no stable solutions. While in this sys- 
tem, the spin-orbit interaction has been considered to 
be relevant in the small or no JT distortion so far, the 
spin-orbital superexchange interactions due to the elec- 
tron transfers turn out to dominate over the spin-orbit 
interaction. Moreover, if we would take the dominance 
of the spin-orbit interaction, there would be a discrep- 
ancy between the calculated energy-difference and that 
estimated from T^. Thus, we conclude that the spin- 
orbit interaction is irrelevant to the origin of AFM(G) 
state in LaTi0 3 , and the experimentally observed re- 
duction of the moment can be attributed to the strong 
itinerant fluctuations in LaTiC>3 instead of the spin-orbit 
interaction. Indeed, a recent neutron-scattering experi- 
ment reveals the spin- wave spectrum of LaTi03 well de- 
scribed by a spin- 1/2 isotropic Heisenberg model on the 
cubic lattice and the absence of unquenched orbital an- 
gular momentum. 12 ^ This indicates that the spin-orbit 
interaction is not effective in this system. Our results 
support these experimental results and suggest another 
mechanism for the emergence of the AFM(G) state. We 
expect that effects which are not treated in our model are 
responsible for its origin. Recently, possible D 3c i distor- 
tion of the TiOe octahedron is examined as a candidate 
for the origin and nature of AFM(G) state in LaTi0 3 . 34 ) 
In this scenario, the spin-orbit interaction is dominated 
by the i2g-level splitting due to the D^d crystal field. 

In addition, we have also studied a magnetic phase di- 
agram in the plane of the GdFe03-type and cZ-type JT 
distortions in order to examine how extent the physics 
of AFM-FM phase transition in strong coupling limit 
survives when the JT level-splitting competes with the 
spin-orbit interaction. According to the obtained phase 
diagram, the description of the phase transition obtained 
by the previous strong coupling approach is well estab- 
lished in the wide range of JT distortion even when the 
spin-orbit interaction is taken into consideration. 
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